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ABSTRACT 

An alternative representation of the quasistatic nonradiative rates of Gersten and Nitzan 
[J. Chem. Phys. 1981, 75, 1139] and Ford and Weber [Phys. Rep. 1984, 113, 195] is 
derived for the respective parallel and perpendicular dipole orientations. Given the dis- 
tance d of a dipole from a sphere surface of radius a, the representations comprise four 
elementary analytic functions and a modified multipole series taking into account residual 
multipole contributions. The analytic functions could be arranged hierarchically accord- 
ing to decreasing singularity at the short distance limit d — > 0, ranging from d~ 3 over 
d^ 1 to ln(d/a). The alternative representations exhibit drastically improved convergence 
properties. On keeping mere residual dipole contribution of the modified multipole series, 
the representations agree with the converged rates on at least 99.9% for all distances, 
arbitrary particle sizes and emission wavelengths, and for a broad range of dielectric con- 
stants. The analytic terms of the representations reveal a complex distance dependence 
and could be used to interpolate between the familiar d~ 3 short- distance and d~ 6 long- 
distance behaviors with an unprecedented accuracy. Therefore, the representations could 
be especially useful for the qualitative and quantitative understanding of the distance 
behavior of nonradiative rates of fluorophores and semiconductor quantum dots involving 
nanometal surface energy transfer in the presence of metallic nanoparticles or nanoanten- 
nas. As a byproduct, a complete short-distance asymptotic of the quasistatic nonradiative 
rates is derived. The above results for the nonradiative rates translate straightforwardly 
to the so-called image enhancement factors A, which are of relevance for the surface- 
enhanced Raman scattering. 



*http: / /www. wave-scattering. com 



1 



1 Introduction 



The distance dependence of the nonradiative rate of a fluorophore or a semiconductor 
quantum dot relative to a metallic nanoparticle (MNP), or a nanoantenna, is of crucial 
importance for a number of novel techniques and nanodevices combining biomolecules, 
quantum dots, and MNP [HI21II2HI51E1IZ1IS1IS]. The MNPs can quench fluorescence 
as much as 100 times better than other quenchers of fluorescence, such as DABCYL, and 
open new perspectives in the use of hybrid materials as sensitive probes in fluorescence- 
based detection assays PQ. In the case of single-stranded (ss) DNA probes labeled with 
a thiol at one end and a dye at the other, the ssDNA molecules self-organize into a con- 
strained conformation on the MNP surface and the fluorophore is completely quenched by 
the particle through a distance-dependent process. Upon target (complementary oligo) 
binding, the constrained DNA conformation is opened because of a dramatic increase in 
the DNA rigidity that becomes double-stranded after hybridization. Thereby the fluo- 
rophore is separated from the particle surface, and the fluorescence of the hybrid molecule 
can increase by a factor of as much as several thousand as it binds to a complementary 
ssDNA pp. This structural change generates a fluorescence signal that is highly sensitive 
and specific to the target DNA [H [2], resulting in a molecular beacon that can detect 
minute amounts of oligonucleotide target sequences in a pool of random sequences and 
provide an improved detection of a single mismatch in a competitive hybridization assay 
for DNA mismatch detection [TJ. MNP acceptors can also be employed for probing changes 
in distances for protein interactions on DNA using a molecular ruler approach involving 
nanometal surface energy transfer (NSET) from optically excited organic fluorophores to 
small MNPs (5j E]. Whereas the traditional Forster resonance energy transfer (FRET) 
involving molecular acceptors and donors is efficient for the acceptor-donor separations 
up to 100 A, the use of MNP acceptors could more than double the traditional Forster 
range up to 220 A [6]. As it becomes apparent below, the reason for the extension of 
the traditional Forster range is a relaxed distance dependence relative to FRET, because 
the particle-fluorophore interaction could no longer be described by sole dipole-dipole 
interaction. 

The nonradiative rates are derived from the total Ohmic loss, or the Joule heating, 
which is in general determined as the power P a t, s absorbed inside the particle [TQl [TTJ [T2j 



where the volume integral extends over entire absorbing region. The quantity Q is the 
steady (averaged) inflow of energy Q per unit time and unit volume from the external 
sources which maintain the field, and which in the gauss units is given by 



where fco is the vacuum wave vector and e" (//') is the imaginary part of the dielectric 
function (magnetic permeability) at the observation point. Alternatively, in the special 
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case of an oscillating dipole source pe luJt , one could consider the dipole moving in a local 
field produced by its surface image and calculate the dipole dissipated power by (see eq 
3.1 of ref [13]) 

P abs = ^lm(p* -E md ) (3) 

where E ind e~ tuJt is the electric field induced by the dipole at the dipole position r d , i.e. by 
the current source 

j(r,t) = -iwpe- fart <y(r-r d ) (4) 
The nonradiative rate T nr is then given by the correspondence principle 

T nr = ^ (5) 

For the sake of definiteness, let the sphere center be located at the coordinate origin, the 
dipole position denoted be described by the position vector r d , r d = \r d \ > a (see Figure 
HJ, and the respective e and e be the sphere and host dielectric constants. Although exact 
expressions for the total, radiative, and nonradiative rates for a dipole interacting with 
a spherical particle pm [HI [THJ, [19], including the case of a multi-coated sphere [HI IT5] . 
are known, and the F77 source codes could be freely downloaded [20], the quasi-static 
results of Gersten and Nitzan [121 EI] and Ford and Weber [13J (cf eqs [6] and [7] below) 
have nevertheless still retained their value. The quasi-static results for small enough sphere 
radius a agree rather well with exact calculations [HI [161 E] an d experiment [7] , and allow 
one for a valuable qualitative analytical insight. For a small separation d of a fluorophore 
from a sphere surface, such that the dimensionless distance parameter 8 = d/a «C 1, the 
quasi-static results predict a d~ 3 dependence of the nonradiative rates (see eqs [32] and [33] 
below). The d~ 3 short-distance dependence corresponds exactly to the case of a dipole 
located at the distance d above a half-space characterized by the dielectric constant e 
[HI [HI [2T] . However, as shown in Figures HI [6] the d~ 3 short- distance dependence begins 
to deviate rapidly from the converged rates for 5 > 0.1 with increasing 5. Thus, for a small 
MNP with a < 10 nm, the d~ 3 short-distance dependence of the quasi-static nonradiative 
rates turns out to be valid only in an immediate proximity to the MNP (d < 1 nm) where 
the very use of the quasi-static rates is questionable. In the present work we remedy 
the above shortcomings by providing an analytic formula, which describes the distance 
behavior of the nonradiative rates over all distances with an unprecedented accuracy. 

In what follows, we derive in Section [2] an alternative superconvergent representation 
of the quasistatic nonradiative rates (eqs HH EH [25] below) with drastically improved 
convergence properties. The superconvergent representation will be obtained by dividing 
the total contribution of each given multipole into a number of partial contributions (eqs 
[T6l [221 below) so that the essential part of the partial contributions could be summed up 
analytically (eqs [T71 [231 below) leaving behind a residual infinite multipole series. The first 
four exact sums are expressed in terms of elementary analytic functions of u — (a/r d ) 2 , 
with each subsequent function having a lesser singularity at the short distance limit d — y 0, 
ranging from d~ 3 over d" 1 to ln(d/a). Convergence properties of the alternative represen- 
tations are examined in Section El On keeping mere residual dipole term of the modified 
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multipole series, the alternative representations are demonstrated to agree with the con- 
verged quasistatic nonradiative rates (eqs El and [7j) of Gersten and Nitzan [121 El] and 
Ford and Weber [13] on at least 99.95% over entire length interval, for arbitrary particle 
sizes and emission wavelengths, and for a broad range of dielectric constants. In Section 
S] the contribution of the analytic terms relative to that of the residual multipole series 
is analyzed. The alternative representations make it also possible to derive a complete 
short distance asymptotic of the classic quasi-static nonradiative rates including all sin- 
gular terms in the distance parameter 5. This is performed in Section El A two-term 
asymptotic will be shown to provide a significant improvement over the conventional d~ 3 
asymptotic over extended length scale. Keeping more terms in the asymptotic expansion 
improves the asymptotic in an immediate sphere proximity 5 <C 0.1, but it comes at the 
expense of worsening the precision for 5 > 0.25. In Section El common approaches for 
fitting experimental nonradiative rates are critically examined in light of the results of the 
present paper. In Section [7J the results for the nonradiative rates are shown to straight- 
forwardly translate to the so-called image enhancement factors A, which are of relevance 
for the surface-enhanced Raman scattering (SERS). The alternative superconvergent rep- 
resentations provide a significant improvement over the familiar d~ 3 short- distance and 
d~ & long-distance behaviors, and could be especially useful for the qualitative and quan- 
titative understanding of the distance behavior of nonradiative rates of fluorophores and 
semiconductor quantum dots involving nanometal surface energy transfer in the pres- 
ence of metallic nanoparticles or nanoantennas [U El El El E] ■ The remaining points are 
discussed in Section El We then conclude with Section El 



2 Superconvergent representations of the quasistatic 
nonradiative rates 

When calculating decay rates for a dipole interacting with a spherical particle, one dis- 
tinguishes the cases of a perpendicular (_L) and parallel (||) dipole orientation relative to 
the sphere surface. The case of a general dipole orientation is then reduced to a linear 
combination of the two particular cases. Let e = s/sq be the relative dielectric constant, 
and let /j, denote the molecular dipole. The nonradiative rates of Gersten and Nitzan 
[T2"l |2"T] for a perpendicular dipole orientation are (cf eq B.24' of ref [2Tj ) 

= -l^y(i+i) 2 - ^^itt (6) 

2he a 3 ^ ' \r d J e + l -f K J 

whereas the rates for a parallel dipole orientation are (see eq B.45' of ref 

I 1 2 00 / \ 2Z+4 1 

MID = -^EIWIPW)^ > m?T¥ 
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yVj + l)(-) Im^^ (7) 

where Im denotes the imaginary part. The respective 2nd lines of eqs [6] and [7] have been 
obtained on using 

Im 1 = ' Im ^l (8) 
which follows from the identity: 

e-l = e + - — (9) 

Occasionally, the Gersten and Nitzan [T2| |2T] expressions are written with an intrinsic 
molecular dipole moment /x that is related to /x by (cf. eqs B.13, B.15 of ref [T2] ) 

/x = jx + aE loc = (10) 

Here Ei oc is the local electric field that is in general affected by the presence of a sphere, 
a is the molecule polarizability, and A is a corresponding image enhancement factor. 
The intrinsic molecular dipole moment /j, is what would be found for a totally isolated 
molecule far away the sphere. A is significant only very close to the surface, in a region 
where the use of a classical approach and of the point dipole model for the molecule are 
questionable. Hence A may be disregarded in computing the actual decay rates [21] . The 
formulas (eqs [6] and [TJ were independently confirmed by Ford and Weber, who determined 
the power dissipated by an arbitrary oriented dipole as (eq 3.47 of ref fl~3~] ) 



abs 



£(/+i) (i+iM+L* («) im y - £ y (ii) 



2e a3 j=j 

where ^x ± and fxy are the respective perpendicular and parallel components of /x = 

When substituted in the correspondence principle (eqlSJ), the result for P a b s leads to the 

nonradiative rates that are identical to those obtained by Gersten and Nitzan [T2l |2T] . 

In what follows, we employ a simple recipe to arrive at an alternative representation 
of the Gersten and Nitzan (GN) expressions (eqs M and [7J) with drastically improved 
convergence properties. First note that the 'Im' sign in eqs [6] and [TJ can be brought in 
front of the summation sign. Afterward, the task of deriving an alternative representation 
of the power series in eq[6]for a perpendicular dipole orientation reduces effectively to that 
for the power series 



S ± = ±(l + D 2 ^t (~ Y = ^ £ (12) 
£r £ + 1 + 7 W e + l^/ + w 



where w and w are the shorthands for 
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In terms of S i 



m 2 T q 

2h^ U lmSl 



(14) 



Now follows a point of crucial importance. This is an iterative decomposition of the coeffi- 
cients of the power series in eq [12] into decreasing powers of I. The iterative decomposition 
could be illustrated as follows. We begin to single out the leading I 2 dependence as follows 



/(Z + l) 2 (l + v-v)(l + lf 



I + V 



I + V 



(l + l 



,2 v{i + iy 



l + V 



(15) 



Note that the last term in eq[15]is effectively of the order I for I 3> 1, which is smaller 
than the I 2 order of the initial term. One can thus repeat the preceding step on the last 
term in eq[15l etc. By continuing the above iterative decomposition, one arrives at 



1(1 + 1) 
l + v 



(l + l) 2 -v(l + l)-v(l-v)- 



v(l-v) 2 v 2 (l-v) 2 v 3 (l-vf 



+ 



l 2 + l(2-v) + (l 
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I I 2 l 2 (l + v) 

v(l-v) 2 v 2 (l-v) 2 v 3 (l-v) 2 



I 



+ 



p 



l 2 (l + v) 



(16) 



On applying the above decomposition in eq[12j the total contribution of each given multi- 
pole is divided into a number of partial contributions. The advantage of the decomposition 
is in that the first four partial contributions in eq. [16] could easily be summed up analyt- 
ically. Indeed, upon substituting eq [16] back into eq[12] one obtains on using elementary 
summation formulas (eqs 146] - 149] of Appendix [AJ 



u 2 S A 



e-1 



e + 1 



u 



u(l + u) u(2 - v) 



1 



u 



1 



u 



u(l 



u 



+ v(l - v) 2 ln(l -u) +v 2 (l -v) 2 F(u,v) 



where the remainder 



F(u,v) = J2 



u 



OO ; 
U 



i in + v) 



u 



i=i 



^l 2 (l + v) 



;is) 



takes into account a residual contribution of the original multipoles that has not been 
comprised by the analytic functions in the square bracket of eqdT] On substituting eqlTTl 
in eqfT41 



r„ r (.L) 



e e 



2heoa 3 

2 



Im 



e + e Q 



u 3 (l + u) 



u 



2e + e 



+ 



[1 — u) 3 (1 — u) 2 e + e 1 — u (e + e 



2^2 



u ln(l - u) + 



u 2 F(u, v) 



(e + e ) 3 ' (e + e )* 

where we have substituted for v in eq [13 according to eq [13] 



(19) 
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The task of arriving at an alternative representation of the power series in eq [7] for a 
parallel dipole orientation reduces effectively to that for the power series 



5,| = 5^ 1(1 + 1) 



1 



i=i 



e+l + j \r d 



21 



l^/ 2 (/ + l) , 



e + 



i=i 



l + v 



-u 



(20) 



where u and v are the same as in eq O In terms of S\\ 

12 



Aheoa 3 



u Im S\\ 



(21) 



When the iterative decomposition described above in connection to S± is applied to the 
coefficients of the power series for Su in eq[20], one finds 



l 2 (l + l) 
I + v 



I 2 + 1(1 - v) - v(l - v) + 



v 2 (l—v) v 3 (l — v) v 4 (l—v) 



(22) 



Upon substituting eq [22] back into eq 
(eqs |4~6H49]) , one arrives at 



I I 2 l 2 (l + v) 

and on using elementary summation formulas 



u 2 5,i 



■ u 



e+ 1 
uv(l — v) 



u(l + u) u(l — V 



1- u)' A ' (1 - u) 

2/i „.\ „.\ „.3/ 



1 



v 2 (l - u) ln(l -u)- v A (l - u) F(u, v) 



(23) 



with the remainder F(u, v) has been defined by eq[T51 Compared to eqUTJ the w-dependent 
terms in eq[23]are, up to their prefactors, identical. On substituting eq|23]m eq| 
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le + go J 



u 2 ln(l 



(g + g ) z 



u F(u,v) 



1 — U(e + Eq) 2 

(24) 



where we have substituted for v in eq [23] according to eq [ 

The respective formulas (eqs [T9l and [241) provide the sought alternative representations 
of the quasistatic nonradiative rates Gersten and Nitzan (T2J [21] and Ford and Weber [TJ] 
for the parallel and perpendicular dipole orientations. No approximation has been used yet 
and the formulas of eqs [19] and [2U are fully equivalent to the original GN expressions (eqs 
[6] and [7]). The respective alternative representations comprise four elementary analytic 
functions together with a residual multipole series (cf eqs [131 E 



u 2 F(u,v) = J2 



g + go 



.1+2 



j^l[le + (l + l)e ] 



u 



(25) 



One can verify that the dependence of the nonradiative rates on e and Eq remains es- 
sentially only through the relative dielectric contrast e, whereas Eq merely multiplies a 
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universal and particle size independent part expressed in terms of the relative parameters 
u and e. The elementary analytic functions in the square brackets in eqs [19] and [2H are 
ordered hierarchically from left to right according to decreasing singularity in the short 
distance limit d — > 0. Indeed, on making the use of dimensionless distance parameter 
5 = d/a <C 1 and on expanding the binomial u = [1 + (d/a)]~ 2 = (1 + 5)~ 2 , one finds in 
the limit 5—^0 

(1 -u)- 3 ~ (25)- j (26) 

The 4th term exhibits a logarithmic singularity ~ In 25. The power series u 2 F has a finite 
limit for 5 — > 0, or equivalently u — > 1: the first series in the second equality in eq [TBI 
sums up to dilogarithm L^w) (see eq 150]) . which is a regular function in the limit u — > 1 
(see eq [52]) , and the coefficients of the second series in the second equality in eq [18] are for 
/ 1 bounded from above by those of Li 2 (w) (cf eq l50]) . 

The reason why the respective iterative decompositions (eqs [TBI and 122]) of the multipole 
contributions in the power series in eqs [T2] and [2U] enable to disentangle various singular 
terms can easily be understood. The power series in eqs [12] and [20] have the radius of 
convergence \u\ = 1, corresponding to d = (cf eq[13]). With the decomposition in eqs [TBI 
or [22] being effectively an expansion of the multipole contributions in eqs [12] and [20] into 
decreasing powers of I, each subsequent term in either eqdBJor eq [22] corresponds to the 
series contributing to a decreasing order of singularity in the limit 5 — » 0. Obviously, the 
described iterative decompositions could, in principle, be continued further also for the 
coefficients of the remainder F(u,v) in eq[18] The alternative representations (eqs [191 and 
121]) of the nonradiative rates will be shown to 

• possess a drastically improved convergence properties; 

• enable a compact analytic description of the GN nonradiative rates over entire pa- 
rameter range; 

• enable one to derive a complete short- distance asymptotic behavior of the GN non- 
radiative rates in the dimensionless distance parameter 5; 

• translate straightforwardly into corresponding results for the image enhancement 
factors A. 

The properties will be dealt with in the subsequent sections. 

3 Convergence properties 

For practical calculations a cut-off l max has to be imposed on infinite power series in each 
of eqs [6] [7] and [25] In what follows, the average nonradiative rates will be plotted. The 
latter are obtained by angular averaging over different dipole orientations at a given dipole 
radial position, resulting in (see eq 7 of ref [T9"] ) 

r = [2T(||)+r(±)]/3 (27) 
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(Note in passing that the rates r(||) enter with a factor 2 into the orientational average, 
because there are two linearly independent parallel dipole orientations possible.) The 
reason behind plotting the average nonradiative rates is that they encompass the particular 
cases of two relative dipole orientations considered in the previous section. The average 
nonradiative rates obtained for a given cut-off l max will be plotted as the ratios to the 
converged nonradiative rates. The converged rates were here and below calculated with 
the cut-off l max = 500 in eqs [6] and in order to achieve a sufficient convergence at an 
immediate sphere proximity. Fortran F77 computer code is freely available on-line [22] . 
It is well known a sphere may exhibit a series of the /-polar localized surface plasmon 
resonances (LSPR), with the respective resonance conditions given by the formula [23J 

e = --^e (28) 

Two particular cases are considered: (i) an off-resonance case, where the real part e' of e is 
well outside the resonance interval — 2e < e' < — e , and (ii) an on-resonance case, where 
e' is close or within the resonance interval while the imaginary part e" of e is relatively 
small. For the sake of simplicity, the host medium will be in what follows assumed to be 
air (eq = 1). 



3.1 Off-resonance case 

Figure |2] displays the so-called average nonradiative rates for the off-resonance case corre- 
sponding to a spherical nanoparticle with the radius a = 10 nm, the emission wavelength 
of the dipole emitter 612 nm, and the dielectric constant at the emission wavelength 
e(612 nm) = —15.04 + 1.02z [T3| ITT] . The dependence of the nonradiative rates on the 
emission wavelength of a dipole emitter is only implicit through e and Eq. Additionally, 
the dependence of the nonradiative rates on the particle size only enters through the pref- 
actor a -3 , which multiplies a universal and particle size independent part expressed in 
terms of the dimensionless parameter u. Therefore, although the above parameters were 
tailored to a silver nanoparticle (AgNP), the behavior shown in Figure H] is expected to 
hold (at least qualitatively) for arbitrary particle sizes and sphere materials, as long as 
one deals with an off-resonance case. Typically, the situation in Figure [2] corresponds to 
a rather common case of a ffuorophore or a quantum dot with the emission wavelength 
of > 480 nm in the case of a AgNP or with the emission wavelength of > 550 nm in the 
case of a gold nanoparticle (AuNP). 

The convergence of the conventional GN representations (eqs [6] and [7]) of nonradiative 
rates is illustrated in Figure [2^. The conventional representations are characterized by 
a rapid drop in precision at an immediate sphere proximity. Even with l max = 20, the 
convergence for d/a < 0.1 is less than 77%. Obviously, the nonradiative rates cannot be 
approximated by taking into account solely the dipole-dipole interaction [TUl HZ]. Com- 
pared to that, the alternative representations of nonradiative rates shown in Figure [2b 
exhibit a superconvergence - already for l max = 1, i.e. by taking into account a mere resid- 
ual dipole contribution in eq |25l the alternative representations agree with the converged 
rates up to at least 99.95% over entire distance interval. 
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3.2 On-resonance case 



An on-resonance case is illustrated in Figure [31 which shows the ratios of the approximate 
and converged average nonradiative rates for a = 100, 6q = 1, the emission wavelength of 
the dipole emitter 340 nm, and £(340 nm) = — 1 + 0.6i. Similarly to Figure [2] the latter 
value of e corresponds to the dielectric constant of AgNP at the wavelength slightly below 
340 nm. However, the situation shown in Figure |3] is not generic, but rather specific to 
a metal with very low imaginary part e" at the wavelength where e' = —1. For instance, 
in the case of Au one would find e(400 nm) = —1.07 + 6.486i [23], i.e. ten-times larger 
e" than in the case of silver, which would rather correspond to an off-resonance case (see 
also Section |3]). 

Compared to the off- resonance case shown in Figure [2^, the dipole contribution shown 
in Figure is significantly reduced at the corresponding values of 5: from « 50% down 
to ~ 30% for 8—1, and from w 80% down to ~ 60% for 5 = 2. Compared to the 
off-resonance case of Figure |2b, the precision of the alternative representations shown 
in Figure |3b has for l max = 1 again a minimum at S « 1. At 5 ~ 1 the precision 
drops from 99.95% in the off-resonance case shown in Figure |2b down to 87% in the on- 
resonance case. The precision is still very good. Note in passing that one customarily uses 
the respective d~ 3 asymptotic behavior (see Figure 13]) and d~ 6 asymptotic behavior (see 
the lmax = 1 contributions in Figures [2^ and |3^) to fit experimental data even if much 
lower precision (less than 80%) is guaranteed. Interestingly, the inclusion of the residual 
quadrupole contributions for l max = 2 rapidly restores the convergence of the alternative 
representations, which becomes not worse than 98% over entire length interval. Therefore, 
the only difference compared to an off-resonance case is that one should increase l max from 
1 to 2 in order to obtain a fairly accurate description of the nonradiative rates. 

3.3 Summary of convergence properties 

It should be not surprising that the alternative representations have better convergence 
properties. After all they were derived by summing exactly the first four partial contri- 
butions of each given multipole in eqs [TH] and |22J up to an arbitrary high order. Since 
higher order multipoles in eqs [6] and [7] become increasingly relevant with decreasing d, 
it was to be expected that the alternative representations would perform much better in 
an immediate sphere proximity relative to eqs |6] and [7] What is surprizing here is an 
astonishingly fast convergence. 

In order to understand the superconvergence, note the following. The coefficients of 
the power series (eq l25]) . which takes into account residual multipole contributions, are of 
the order 0(l~ 2 ), whereas the coefficients of the original multipole series (eqs E] and [7]) are 
of the order 0(l 2 ) for I ^> 1. Therefore, beginning with the residual quadrupole (I = 2) 
contribution, the residual contribution of higher order multipoles in the series (eq [25|) has 
been reduced by the factor of l~ A compared to the original GN expressions (eqs [6] and 
[7]). Of course, the missing part of the multipole contributions did not disappear - it has 
become comprised in the analytic terms. The effect of higher-order multipoles comprised 
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in the contribution of the elementary analytic functions in the square brackets in eqs [19] 
and [2H is appreciable over entire distance range. Indeed, the alternative representations 
shown in Figure [2b agree with the converged rates within 0.03% for 5 ~ 2 and l max = 1, 
whereas a mere ~ 80% precision is achieved by the dipolar term in the conventional 
GN representation (eqs E] and [7]) of the nonradiative rates shown in Figure [2h- Thus 
the drastically improved convergence of the alternative representations (eqs dH] [2H, |2"51) 
can be viewed as a consequence of the fact that the elementary analytic functions of the 
alternative representations (eqs [19[ l24"j) comprise an essential part of the contribution of 
an infinite number of higher order multipoles. A slightly worse convergence properties in 
an on-resonance case is then caused by an increase of the quadrupole contribution relative 
to the dipole one. The reason of that increase will be further explained in Section H] 

4 The contribution of the analytic terms relative to 
that of the residual multipole series 

In turns out expedient to determine of how much of the higher-order multipole contribu- 
tion is comprised in each of the analytic terms in the square bracket of eqs [19] and [211 For 
that purpose we have considered various approximations to the average nonradiative rates, 
which were obtained by keeping a gradually increasing number of the analytic terms. The 
results for an off-resonance case corresponding to Figure [2] are plotted in Figure HI It is 
obvious from Figure H] that keeping mere first three terms in the square brackets in eqs 
[19] and [24] reproduces the converged rates already within 10% over entire distance range, 
whereas keeping the first four terms in the square brackets in eqs [19] and EH reproduces 
the converged rates within 1.2% over entire distance range. We recall that the last but 
one series in eq [TBI sums to dilogarithm Li 2 (w). One can then substitute 



into each of eqs [Tj2 and EU and consider an approximation which results by neglecting the 
last series, i.e. by approximating F(u,v) by Li2(w). Such an approximation [comprising 
all four singular terms together with Li2(w)] then reproduces the exact results within 
remarkable 0.12% over entire distance range. A similar agreement has also be obtained 
for a = 0.5 nm AgNP as shown in Figure 

The results shown in Figures @] and |5] suggest that our analytic terms have to comprise 
a substantial part of the dipole contribution. Indeed, the dipole {l m ax = 1) contribution 
accounts for 50% and ~ 80% of the average nonradiative rate for 5 = 1 and 5 = 2, 
respectively, in Figure Ek- Therefore, without taking into account an essential part of 
the dipole contribution, less than 50% agreement with the nonradiative rates would be 
possible for 5 ~ 1 in an off-resonance case. By keeping the complete dipole contribution 
one could always achieve a good agreement with the exact rates for sufficiently large d. In 
the latter case the sphere response is that of a polarizable point, and the long-distance d~ 6 




(29) 
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dependence is that corresponding to a dipole-dipole interaction [13]. Indeed, for r d ^> a, 
so that one can approximate with d, the distance dependence of the nonradiative decay 
is obviously dominated by the I = 1 dipole term in the GN series (eqs [6] and [7J , and the 
nonradiative decays follow a d~ 6 distance dependence. Interestingly, the square brackets 



in eqs [19] and [2H exhibit a leading d 



~ u 



behavior in the long distance limit, in 



which case m C 1. Therefore, as long as the residual dipole contribution in the remainder 
F(u,v) (defined by eq 1251) in eqs fT9l and 1241 is negligible compared to the original dipole 
contribution in eqs [6] and [7J our analytical terms are expected to provide very good 
approximation to the nonradiative rates. Note in passing that the product u 2 F(u,v) in 



eqs [19] and [2U reduces for l Ti 



1 to (cf eq[25 



u F(u,v) -» 



e + e 



u 3 (30) 

Thus the condition that the residual dipole contribution in F(u,v) is negligible compared 
to the original dipole contribution in eqs M and [7J translates into 



Im 



2e r 



Im 



e-e 



e + 2e 



where c 



2^2 



for a perpendicular dipole orientation and c 



(31) 



£q£ for a parallel dipole 



orientation. Obviously, if the values of e and e are spaced apart (e.g. e = —15.04 + 1.02i 
and Eq = 1), the left-hand side of eq [31] is strongly damped by the denominator of the 
first fraction (e.g. by w 14 2 = 196 in the examples shown in Figures [2] H] and [5]), and the 
condition (eq IST]) is satisfied. 

The condition in eqEHenables one also to determine a worst case scenario, in which case 
the original multipole contributions comprised in the analytic terms of our superconvergent 
representations (eqs [T9l and |24"|) are at a minimum. The latter would occur for e m —So, 
in which case the denominator of the first fraction on the left-hand side of eq [31] would 
work against the inequality [3U Because the residual multipole contributions in the series 
(eql2"5T) have been reduced by the factor of Z -4 compared to eqs [6] and [7J a violation of the 
condition in eq[3T]would not affect much high order multipoles. Typically, only the relative 
residual contribution of a few low order multipoles will be affected. This is exactly the 
scenario that has been illustrated in Figure [3] where it has been shown necessary to keep 
the residual quadrupole contribution in the remainder F in order to ensure a reasonable 
convergence. 



5 Short-distance asymptotic behavior 

With decreasing d, higher order multipoles in eqs O and [7J become increasingly relevant. 
Ultimately, in the limit of a sufficiently small dipole-sphere separation 5 = d/a < 1, an 
infinite number of multipoles in eqs [6] and [7J contributes to the change of the d~ 6 long 
distance dependence into the d~ 3 short distance dependence (see eq B.24" of ref [2T] ) 

r nr (i_) -> Im = Im £ -^- (32) 

v ' Ahe Q d? e + 1 8he d 3 e+1 K ' 
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and (see eq B.45" of ref [2TJ) 



r «r(ID ^ - q^s Im -TT = T7^ Im "-TT C 33 ) 
8ne d 6 e + 1 16n,£:o« £ + 1 

(The respective 2nd equalities in eqs [32] and [33] have been obtained on using eq M in the 
limit I — y oo.) For a comparison, the power dissipated by an oscillating dipole at the 
distance d above a planar interface of a metallic half-space is (cf eq 3.23 of ref 



abs 



^_ L? + I M 2 \ Im £^£o (34) 



One can easily verify that upon substituting P a b s into the correspondence principle (eq 
EJ, the limit values of nonradiative rates given by eqs [32] and [33] are recovered. Thus the 
leading d~ 3 short- distance dependence corresponds exactly to the case of a dipole located 
at the distance d above a half-space characterized by the dielectric constant e. 

Given the alternative representations (eqs [19] and El]), it is possible to derive the 
short-distance asymptotic behavior of the Gersten-Nitzan expressions involving all singu- 
lar terms in the limit 8 — >■ 0. A slight nuisance is that the first term in the square bracket 
in eqs [T71 and |2"31 contributes also to the 5~ 2 and 5~ l terms (see eq[53]). Similarly, the sec- 
ond term in the square bracket in eqs [T7] and [23] contributes an additional 5 _1 term (see 
eq l54l) . After taking into account the sub-leading singular terms according to elementary 
formulas (eqs 15311551 of Appendix [A} and substituting back into eq [TTJ one finds in the 
limit 5 -» 

u 2 S ± ~yL + ^ + ^± + aiog \n{25) + 0(1) (35) 



where 



0-3 



a_ 2 



a_i 



e-eo 

4(e + £ )' 
(e - e )(e + 3e ) 
8(e + e ) 2 
(e-e )(e 2 + 3el) 
8(e + e y 



e e 2 (e-e ) 

aio 9 = (£ + £o)4 (36) 
On repeating the steps leading from eq [IT] to eq [35], one finds in the limit 5 — > 

u2S \\~ b -^ + b -^ + b -f + 6 W k(25) + 0(1), (37) 

where b_ 3 = a_ 3 and 

(e - e )(3e + 5e ) 
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'log 



(e - s )(3e 2 + 8e e + 9eg) 

(38) 



8(e + e )z 
e%e(e - e„) 



(e + 



The corresponding short-distance asymptotic of the nonradiative rates follows then straight- 
forwardly on substituting eqs [35] and [37] into eqs [14] and EH respectively. As a consistency 
check, the term proportional to <5~ 3 could be shown to reproduce eqs [32] and [33] respec- 
tively. Looking at the coefficients (cf eqs [36] and [38]) of the asymptotic expansions, it is 
obvious that by a judicious choice of the dielectric constants e and e one could switch off 
either 5~ 2 or S^ 1 terms in eqs I3"5l and I3T1 For instance, 

=>• a__2 = 

==>- a_i = 

6_i = (39) 

A short-distance asymptotic behavior of the average nonradiative rates is shown in Figure 
[6] The configuration is the same as in Figure [3] Unfortunately, the use of the above 
asymptotic expansion is limited to a rather short distance interval. Indeed, the precision 
drops to below 90% and 80% for 8 ~ 0.15 and 8 ~ 0.2, respectively. The major effect 
limiting the applicable distance range is a limited range 8 < 0.2 of the validity of the 
asymptotic expansions (eqs 1531155]) of the elementary analytic functions in the square 
bracket in eqs [19] and [24] [25]. Note in passing that the two-term asymptotic shows overall 
the best properties for an improved fitting, because it provides a significant improvement 
over the conventional d~ 3 asymptotic over extended length scale. Keeping more terms 
in the asymptotic expansions (eqs [35] and l3Tj) improves the asymptotic in an immediate 
sphere proximity 8 <C 0.1, but it comes at the expense of worsening the precision for 
8 > 0.25. 




6 A comparison with the plane-surface d 3 distance 
dependence 

In a large body of current literature one attempts to fit experimental nonradiative rates 
with a power law dependence l/d a , where a — 3 for 8 < 1, and a is between 3 and 4 when 
the separation satisfies l<5<4[2][5][6l[8][9l [26] . However, a d~ 3 distance dependence 
(eqs 132] and 133]) is typically limited to much smaller distance d/a < 0.1 (see Figured]) 
and its precision rapidly decays with the distance. Indeed, the leading d~ 3 asymptotic 
equals twice the converged rates for 8 ~ 0.5 and four-times the converged rates for 8 ~ 1 
(not shown) for the set-up considered in Figure [6] If one wanted to remain within 10% 
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of the GN result, one is limited to 5 < 0.1, which implies sub-nanometer distances for 
MNP radii a < 10 nm. This is clearly not only unsatisfactory to reliably describe recent 
experiments [H [21 El El [6] , but also the very use of the quasi-static rates is questionable 
for such small distances. It is obvious from Figure El that a two term asymptotics yields a 
substantial (essentially 100%) improvement in precision over the leading d~ 3 short distance 
dependence: one can remain within 10% of the GN result for 5 < 0.18 and within 50% 
of the GN result for 5 < 0.5. The two-term asymptotic then equals twice the converged 
rates for 5 ~ 1 (not shown). Another example of a substantial improvement of a leading 
asymptotic by a subleading term has been discussed in ref (27] (see Figure 4 therein). 

In an intermediary region of distances 1 < 5 < 4 there appears to be no theoretical 
basis for any pure power law dependence such as l/d a , where a is between 3 and 4 
[21 El El El E]- Obviously, as shown in Figure El even a mixture of different negative powers 
of 5 in full asymptotic expansions (eqs. I3"5l and 157]) yields rather poor approximation to the 
nonradiative rates. A significant improvement in describing the distance dependence of the 
nonradiative rates can instead be achieved by using the singular analytic terms (see Figures 
HIE]). In view of its simplicity, unmatched precision, and easy use it is therefore preferable 
to use our superconvergent representations (eqsEHl EU E5]) with a cut-off l max = 1, or in 
some cases l max = 2, as fitting formulas for the experimental nonradiative rates. 

7 Image enhancement factors 

According to eq [TOj the value of A describes the change in the net molecular dipole 
moment /x, which in turn determines the relevant scattering cross sections. Therefore, the 
short distance behavior of the respective image enhancement factors A is of interest for 
the surface-enhanced Raman scattering (SERS). Given the respective definitions of S± 
and S\\ (eqs [T21 and |2"0~|) . the image enhancement factors for the perpendicular and parallel 
dipole orientations are (after a trivial recasting of eq B.16' of ref [2T] ) 

I — 1 

and (see eq B.42' of ref (21]) 

A|| = A V Kl + 1) ( ~) = - 2 $\\ (41) 

Therefore, the results obtained in preceding sections for the nonradiative rates translate 
straightforwardly to those for the image enhancement factors by substituting the relevant 
expressions for u 2 S± and u 2 S\\ into eqs E0] and EH respectively. 

8 Outlook 

Our various approximation were extensively compared against the converged quasi-static 
GN rates (eqs El and [7J of Gersten and Nitzan [121 EI] and Ford and Weber [T3]. The 
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approximation could therefore be used whenever the Gersten and Nitzan theory 
applies. In some case one would probably be required to take into account the effect 
of size corrections to the bulk dielectric function [17], which can be straightforwardly 
incorporated. 

We have only considered the nonradiative rates of fluorophores and quantum dots rel- 
ative to a MNP. However, recently a novel class of composite superparticles have been 
introduced, which possess so-called localized dielectric resonances (LDR) having a qual- 
ity factor comparable to that of LSPR [28]. The superparticles open a new avenue of 
applications for the radiative and nonradiative decay engineering, which will be followed 
elsewhere. 

On using a relation between polarizabilities and scattering T-matrix elements, one 
could contemplate to include a dynamic depolarization and radiative reaction terms [29] 
and consider a generalization of multipolar polarization factors according to 



l-(e + l) 



2(21+3) 



lJ r- r\ 1+1 I x2 ( m 3(21+1) - | Q+l)(2Z+3) \ ■ (l+l)x*+i /- -p 
' fc " r I " r 21(21+3) \ lfc 21-1 fc "r" 21-1 I 'i(2i-l)!!(2Z+l)H V fc L ' 



(42) 



where x = koa is the conventional size parameter. Such a generalization is irrelevant 
for small particles [17] . but it would be of interest for larger particles (a > 80 nm), 
where depolarization (~ x 2 term in eq l4"2l) and radiative reaction become increasingly 
important [29]. However, to implement it in all orders, in order to obtain a generalization 
of our super convergent representations, may be difficult. Nevertheless, it turns that an 
analogous generalization of mere dipole order could already provide a sizable improvement 
for larger particles (a > 80 nm) [16] and may sufficiently agree with exact electrodynamic 
calculations pUl dU HHJ HHJ [19] . In this regard note that the results of Merten et al [32] 
could be improved by using a proper dynamic depolarization term that yields a correct 
dipole LSPR position up to the x 2 -order |29j. 

The nonradiative rates are often plotted normalized with respect to the radiative rates 
of a dipole in an infinite host medium. On substituting the Larmor formula for the total 
dipole radiative power 



r (43) 



3c 3 

into the correspondence principle (eqEJ), the dipole radiative rate in the host medium in 
the sphere absence is 



r . (44) 



_ ^ 3 |m| 2 
r 3hc* 

Therefore, the normalization has the net effect of replacing the prefactor in the absolute 
nonradiative rates (eqs El [7J [TH [24]) according to 



heoa 3 x 3 
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9 Conclusions 



In a large body of current literature one attempts to fit experimental nonradiative rates 
with a power law dependence 1/cP, where a — 3 for 5 < 1, and a is between 3 and 4 when 
the separation satisfies 1 < 5 < 4 [21 [5j [6j El [9] . We have shown that such an approxima- 
tion could be highly imprecise (see Figures HI [6b). Instead, alternative super convergent 
representations (eqs[19l |2U [25]) of the quasistatic nonradiative rates of Gersten and Nitzan 
[T2"j |2"T] and Ford and Weber [T3] were derived, which could be used for highly precise, 
simple, and efficient analytic description of the rates. Given the distance d of a dipole 
from a sphere surface of radius a, the representations comprise four elementary analytic 
functions and a modified series (eq [25]) taking into account residual multipole contribu- 
tions. The analytic functions could be arranged hierarchically according to decreasing 
singularity at the short distance limit d — > 0, ranging from d~ 3 over d^ 1 to ln(<i/a). In the 
opposite long distance limit, the analytic functions exhibit a leading d~ 6 behavior. On 
keeping mere residual dipole contribution of the infinite series (eq l25l) . the representations 
typically agree with the converged rates over at least 99.9% over entire length intervals, 
for arbitrary particle sizes and emission wavelengths, and for a broad range of dielectric 
constants (see Figure [2b). The origin of the super convergence was identified and ex- 
plained. The analytic terms of the representations reveal a complex distance dependence 
and could be used to smoothly interpolate between the familiar d~ 3 short- distance and 
d~ 6 long-distance behaviors with an unprecedented accuracy (see Figures H] [5]) . Therefore, 
the representations could be especially useful for the qualitative and quantitative under- 
standing of the distance behavior of nonradiative rates of fluorophores and semiconductor 
quantum dots involving nanometal surface energy transfer in the presence of metallic 
nanoparticles or nanoantennas. As a byproduct, a complete short-distance asymptotic 
of the quasistatic nonradiative rates was derived. The above results for the nonradiative 
rates translate straightforwardly to the so-called image enhancement factors A, which are 
of relevance for the surface-enhanced Raman scattering (SERS). 

The radiative rate in the GN theory is given by a single dipole term (cf eqs B.18' and 
B.43' of ref [21J). The total rate in the GN theory is then obtained by adding the dipole 
radiative rate term to the nonradiative rate. Therefore our results extend straightfowradly 
also to the total rate in the GN theory. 



10 Supporting Information Available: 

You can download the data of the plots and view the plots in different scales by editing op- 



tions of the corresponding Origin^ projects available on-line at http:/ /www. wave-scattering.com/gn. htm 
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A Elementary summation and expansion formulas 

£ „' = ^- (46) 

1=1 

S>'=(T^ (47) 

^ (1_ U )3 

X)y = -Ml-«) (49) 
z=i 4 

and (see eq 5.2.5.4 of ref [30] ) 

> — = - / cfe = Li 2 (u) (50) 

i=i ' Jo s 

defines the so-called dilogarithm function. The identity in eq|48]follows from (see eq 5.2.2.9 
of ref [30]) 

Ec^r^jBw (51) 

i=o K ' 

The identities in eqs H7J and [48] can also be obtained on successively applying the operator 
\u{djdu)\ on both sides of eq[46] The dilogarithm value at the limiting case of u = 1 is 
finite and related to the Riemann zeta function 

Li 2 (l)=C(2) = y (52) 
In the limit 5 — > one finds upon invoking the binomial series 

u = (i + S)- 2 - 1 - 25 + 35 2 -45 3 + 55 4 -G5 5 = l-25{l-q) 

u 2 = (l + <5)- 4 ~ i - 45 + 105 2 - 205 3 + 355 4 

v? = (1 + 5)- 6 ~ 1 - 65 + 215 2 - 565 3 + 1265 4 

u 4 = (1 + ,J)- 8 ~ 1 - 85 + 365 2 - 1205 3 + 3305 4 

u 3 (l + u) ~ 2 - 145 + 575 2 - 1765 3 + 4565 4 

(l-u)- 1 ~ -i (i + ^5 + -S 2 --5 s 
K ' 25 V 2 4 8 

(!_„)-» ~ _L( 1 + 3 * + H^ + ^ 

x 3 1 /, 9. 15. 21 , 3 



and 



u 3 (l +u) 
(1-u) 3 

u 3 
(1-u) 2 
u 3 



1 



1 

A5 3 
1 

452 



3 

45 



23 
16 



+ 0(5 2 



u 



1 9 49 r 

I 5 

25 4 8 



209 
~16 



5 2 + 0(5 3 ) 



(53) 
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Figure captions 



Figure 1 - An illustration of the geometry of the problem. Sphere of radius a is located 
at the coordinate origin. The dipole is positioned outside the sphere at the distance 
r d — l r d| > a from the sphere center. 

Figure 2 - Convergence properties of the conventional GN representations (eqs|6l[7]) of the 
nonradiative rates (a) and novel super convergent representations (eqs [19], EH 125]) of the 
nonradiative rates (b). The average nonradiative rates (eq !27p were calculated for the case 
of a spherical nanoparticle with radius a = 10 nm in air (e = 1), the emission wavelength 
of the dipole emitter 612 nm, and the dielectric constant at the emission wavelength 
e(612 nm) = —15.04 + 1.02i [13] ITT] . The parameters correspond to a silver spherical 
nanoparticle (AgNP), but an analogous behaviour is expected in any off-resonance case. 
The nonradiative rates are plotted as the ratios of the approximate to converged rates. 
The latter were calculated with the cut-off l max = 500 in eqs El [3 in order to achieve a 
sufficient convergence at a sphere proximity. Note different scale on the ordinate axis for 
(a) and (b). 

Figure 3 - Convergence properties of the conventional representations (eqs El E]) of the 
nonradiative rates (a) and novel representations (eqs HHJ), I2"4]) . [23]) of the nonradiative 
rates (b) in an on-resonance case. The average nonradiative rates (eq l2Tj) were calculated 
for the case of a silver spherical nanoparticle (AgNP) with radius a = 100 nm and the 
emission wavelength of the dipole emitter 340 nm. The dielectric constant of AgNP at 
the emission wavelength was taken to be e(340 nm) = — 1 + O.62. Note different scale on 
the ordinate axis for (a) and (b). 

Figure 4 - Comparison of various approximations to the orientationally averaged non- 
radiative rates (eq [27} obtained by keeping a gradually increasing number of analytic 
terms from left to right in the square bracket in eqs []1|] and [24] The average nonradiative 
rates are plotted as the ratios of approximate to converged rates for the case of AgNP 
with radius a = 10 nm in air (sq = 1), the emission wavelength of the dipole emitter 
612 nm, and e(612 nm) = —15.04 + 1.02z (the same configuration as in Figure [2]). The 
plane-surface d~ 3 distance dependence is also shown for a comparison. Keeping the first 
four terms in the square bracket in eqs [19] and [24] reproduces the exact results within 1.2% 
over entire distance range. Upon approximating F(u,v) by the dilogarithm term Li 2 (u) 
(eq |2"9"]). the exact results are reproduced by resulting five analytic terms within 0.12% 
over entire distance range. The last two approximations essentially overlay into a single 
line. 

Figure 5 - The same as in Figure [4] but for a AgNP of radius a = 0.5 nm (e = 1, the 
emission wavelength 612 nm, and e(612 nm) = —15.04 + 1.02z). 

Figure 6 - A short-distance asymptotic behavior of the average nonradiative rates for the 
case of AgNP with radius a = 100 nm in air, the emission wavelength of the dipole emitter 
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340 nm, and e(340 nm) = —1 + 0.6* (the same configuration as in Figure [3]). The Figure 
displays various approximations obtained by keeping an increasing number of terms in eqs 
1351 and [37] The converged GN rates are also shown for a comparison. Obviously plotting 
in the logarithmic scale hides to a large extent the differences with the converged GN 
rates. A comparision of the analytic expressions in eqs [19] and [24] with the inverse cube 
d~ 3 plane-surface distance dependence has been shown in Figure HI 
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